Revisiting the population of Galactic open clusters 
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Abstract. We present results of a study of the galactic open cluster population based on the all-sky catalogue ASCC-2.5 
(I/280A) compiled from Tycho-2, Hipparcos and other catalogues. The sample of optical clusters from ASCC-2.5 is complete 
up to about 850 pc from the Sun. The symmetry plane of the clusters' distribution is determined to be at Zq = -22 ± 4 pc, 
and the scale height of open clusters is only 56 ± 3 pc. The total surface density and volume density in the symmetry plane 
are 2 = 114 kpc~ 2 and D(Zq) = 1015 kpc~ 3 , respectively. We find the total number of open clusters in the Galactic disk 
to be of order of 10 5 at present. Fluctuations in the spatial and velocity distributions are attributed to the existence of four 
open cluster complexes (OCCs) of different ages containing up to a few tens of clusters. Members in an OCC show the same 
kinematic behaviour, and a narrow age spread. We find, that the youngest cluster complex, OCC 1 (log / < 7.9), with 19 deg 
inclination to the Galactic plane, is apparently a signature of Gould's Belt. The most abundant OCC 2 complex has moderate 
age (log / as 8.45). The clusters of the Perseus- Auriga group, having the same age as OCC 2, but different kinematics are seen in 
breaks between Perseus- Auriga clouds. The oldest (log / as 8.85) and sparsest group was identified due to a large motion in the 
Galactic anticentre direction. Formation rate and lifetime of open clusters are found to be 0.23 ± 0.03 kpc~ 2 Myr 1 and 322+31 
Myr, respectively. This implies a total number of cluster generations in the history of the Galaxy between 30 to 40. We estimate 
that less than about 10% of the total Galactic stellar disk population has ever passed an open cluster membership. 

Key words. Galaxy: disk - Galaxy: kinematics and dynamics - open clusters and associations: general - solar neighbourhood 
- Galaxy: stellar content - Galaxy: structure 



. ^ ' 1. Introduction rately than for a single field star (roughly by a factor of yfn 

, where n is the number of cluster members). Actually, they are 

rr* , Open clusters play a double role in astrophysical studies. Since , .. . £1 , . . P , .. ,. , 

?H better tracers of large scale structures of the Galactic disk pop- 
d ' they consist of stars which have been born and lived together, , .. c . . , . . ,. c ., 

. . . i 3 b ulation than held stars. The most comprehensive studies of the 

they provide natural laboratories to prove or rule out theories ^ , .. , , , .. , . OA ,,, T -rrmrsi 

3 F F Galactic cluster population are about 20 years old (Lynga 1982 

of the formation and evolution of stars. Due to the fact that all T . , lr , o. n. u a »l u t a * i 

Janes et al. 1988 1. They were based on the best data avail- 
members have about the same distance from the observers, that , , . . ,. T , ^ , , t ^ , „ . 

able at that time, the Lund Catalogue of Open Cluster Data 

they show common space motion, and have uniform chemical /T - i lr . -, i , „. . , „ . , s , . , . 

1 r (Lynga 1987 hereafter, the Lund Catalogue) and its subset 

composition, each single cluster represents excellent empirical c , . .., n , , . , T „ i lrioo i 

r > t> f r of clusters with 3-colour photometry (Janes & Adler 1982 1, 

reference sequences which are basic for many contemporary , ... . • . . . 

M 1 F 1 Although these studies present an important step in our un- 

scales and calibrations used in astrophysics. ■ . ■• f . u , .. f . u , . , .. 

F 3 derstanding of the general properties of the cluster population, 

On the other hand, studying the open cluster population thgy suffer from incompleteness of the cluster samples ^ 

provides important information that contributes to a better un- from inhomogeneity of the clusters parame ters. 
derstanding of the structure and evolution of the Milky Way. 

In principal, basic parameters like distance, motion, age, and This can be illustrated b Y the following statistics. About 
metallicity can be determined for an open cluster more accu- 1200 clusters were known in the Lund catal ogue by 1988. Only 
400 of them had accurate, but heterogeneous UB V photometry, 

Send offprint requests to: R.-D. Scholz and photometric distances, reddening and age values. Although 



tor almost all clusters apparent diameters were given in the 
Lund catalogue (estimated by eye from sky charts or defined 
by the size of detector's field of view), only about 100 clus- 
ters were studied in a systematic way by use of stellar counts 
(Danilov & Seleznev I1993> . Kinematics, being the most tra- 
ditional field in open cluster studies, was also hampered by 
heterogeneous proper motions and radial velocities (RVs), and 
less than 100 clusters had proper motions reduced to the funda- 
mental reference system (van Schewick ll971> . Space velocities 
were available for a few tens of clusters only. Almost nothing 
was known about characteristic values of cluster statistics, how 
complete and representative the sample was, i.e. the basic pa- 
rameters needed for a systematic study of typical properties of 
the Galactic cluster population. These difficulties were mainly 
due to the fact that data on various clusters were strongly inho- 
mogeneous since they were obtained by different instruments, 
detectors and techniques, with little effort to reduce them into 
one single system. 

Since that time considerable improvement has been 
achieved in the field of new and homogeneous observations 
which could be used for the determination of basic cluster pa- 
rameters and for a systematic search for new clusters. First of 
all, there are the accurate all-sky surveys like Hipparcos/Tycho 
in the optical and 2MASS in the infrared. The first efforts to 
look for open clusters in these data led to the discovery of new 
objects both located in the Solar neighbourhood, and at large 
distances from the Sun. Bica et al. ( 2003a I compiled a list of 
276 embedded infrared star clusters and stellar groups from the 
literature and were able to increase the number by 346 new 
compact and remote objects associated with nebulae by making 
use of the 2MASS data (Dutra et al. l2003l Bica et al. l2003b> . 
Also, nearby clusters were discovered in the Hipparcos/Tycho 
catalogues. Platais et al. ( 1998 ) found 15 nearby, very loose and 
extended clusters and associations in the Hipparcos catalogue, 
whereas Alessi et al. ( 2003 1 discovered 1 1 new clusters within 
0.8 kpc with the Tycho data. 

However, with respect to the astrophysical parameters for 
clusters, the situation had not much improved during recent 
years: the parameter set is neither complete nor homogeneous. 
The on-line list of open cluster data by Dias et al. (2004 
DLAM hereafter) which can be considered as a continuation 
of the Lund Catalogue now contains by a factor of 1.5 more 
clusters than its predecessor, but the degree of completeness of 
this list is still unknown. Since the cluster data in the DLAM 
list are taken from the literature (or from private communica- 
tions), the sets of the derived parameters differ from cluster to 
cluster. Also, the parameters themselves are based on heteroge- 
neous observations and different methods of cluster definition 
and of parameter determination. Whenever using these data for 
cluster population studies, one would meet problems caused 
by uncertain cluster statistics and data heterogeneity. Even one 
of the most uniform lists of cluster parameters, the Loktin et 
al. (|5001 2004 1 catalogue, is based on non-uniform photomet- 
ric data and uncertain membership. In other words, for studies 
of the general properties of the Galactic cluster population, we 
cannot simply benefit from a larger sample of open clusters as 
long as we cannot estimate its homogeneity. 



Ideal preconditions tor cluster population studies are deep 
and uniform sky surveys of sufficiently accurate astrometric 
and photometric data supplemented by radial velocities, metal- 
licities, and spectral classifications. Statistical properties of a 
sample of clusters identified with these data could be well es- 
timated since its completeness should correlate with the com- 
pleteness limit of the surveys used. Applying uniform criteria 
and methods of membership and parameter determination, one 
would get a homogeneous set of basic cluster parameters. In 
the near future an important step can be undertaken in this di- 
rection when the 2MASS data will be combined with the fi- 
nal version of the UCAC catalogue, and with the RAVE and 
SEGUE surveys. In the further future, the observations of the 
GAIA satellite will provide data at an even higher level of ac- 
curacy and homogeneity. 

At the moment, more or less complete and accurate data 
are available only for relatively bright stars and consequently 
for relatively nearby clusters. Nevertheless, it is worth trying to 
use the present data for the re-determination of cluster mem- 
berships and for the re-estimation of basic cluster parameters 
by applying uniform methods to homogeneous date sets. This 
gives us a basis to increase the reliability of the results on 
general properties of clusters, to estimate the possibilities and 
constraints of the data available, and therefore, to update our 
knowledge on the Galactic cluster population, now. 

We started a long trip from the available observations to 
study the population of Galactic open clusters with the com- 
pilation of a complete survey of stars and with the reduction 
(if necessary) of the data to common reference systems. This 
effort resulted in the All-Sky Compiled Catalogue of 2.5 mil- 
lion stars (ASCC-2.5 Kharchenko 2001 ) with absolute proper 
motions in the Hipparcos system, with B, V magnitudes in 
the Johnson photometric system, and supplemented with spec- 
tral types and radial velocities if available. The ASCC-2.5 was 
used to identify known open clusters and compact associations 
from the Lund Catalogue, the Dias et al. (2004i on-line data 
collection, and the Ruprecht et al. J198U list of associations. 
In the ASCC-2.5 we f ound 52 of about 1700 known clus- 
ters (Kharchenko et al. 2004b hereafter Paper I) and discov- 
ered 130 new open clusters (Kharchenko et al. 2005b hereafter 
Paper III). A pipeline was developed to determine cluster mem- 
bership based on kinematic and photometric criteria as well as 
to obtain a uniform set of cluster structural, kinematic and evo- 
lutionary parameters (see paper I and Kharchenko et al. 2005a 
hereafter Paper II). 

In this paper we use the results of Papers II and III to 
study unbiased properties of the cluster population of the lo- 
cal Galactic disk. In §[2] we briefly describe the data set. The 
completeness of the cluster sample and the spatial distribution 
of open clusters in the Galactic plane and perpendicular to it is 
discussed in § The kinematics of the system of open clus- 
ters is considered in § |4] Significant irregularity in the spatial 
and kinematic distributions of clusters revealed complexes and 
groups of clusters which are discussed in §|5] In §[6] we deter- 
mine typical lifetime and formation rate of clusters as it follows 
from the present statistics. Concluding remarks are given in §0 

1 ftp : //cdsarc . u-strasbg . fr/pub/cats/I/280A 



Fig. 1. Distribution of the surface density S of open clusters versus their distance d xy from the Sun projected onto the Galactic 
plane. The dotted lines indicate the completeness limit, the dashed horizontal lines correspond to the average density of "field 
clusters" (i.e., not contained in cluster complexes, s. §0. The bars are Poisson errors derived from cluster counts. Left panel: 
distribution of all clusters. Right panels: distributions of clusters with different ages. The panel designations from Tl to T4 are 
the identifiers of four subsamples of different ages as defined in the text. 



2. Summary of the data used 

The present sample contains 520 already known clusters which 
we identified in the ASCC-2.5 (Paper I) and 130 new open clus- 
ters which we newly detected in the ASCC-2.5 data (Paper III). 
This sample of 650 open clusters is homogeneous in many as- 
pects: it is based on a practically uniform data set (i.e. the cata- 
logue ASCC-2.5 with proper motions and magnitudes reduced 
to uniform astrometric and photometric systems), on a uni- 
form approach of membership selection (kinematic and pho- 
tometric criteria, see Paper I) and on a uniform technique of 
the determination of cluster parameters applied in Papers II 
and III. For each cluster in our sample, we determined dis- 
tance and celestial position (i.e. the location in 3D-space), red- 
dening, size (including core and corona radii), kinematics (the 
mean proper motions in the Hipparcos system), and the cluster 
age computed from ages of individual members near the Main 
Sequence turn-off. On the basis of a new cluster membership 
determination, we derived radial velocities and spatial velocity 
vectors for 359 clusters of our sample. For details of the above 
parameter determination we refer the reader to Paper II. The 
results are presented in two data sets: The Catalogue of Open 
Cluster Data and its Extension 1 . The catalogues are only avail- 
able in electronic form at the CDS online archive 2 . Hereafter, 
we refer to both catalogues with the same identifier COCD. 

The two closest clusters, the Hyades and Collinder 285 (the 
UMa cluster), are missing in our list. Due to their proximity 
to the Sun, they occupy large areas on the sky and require a 
specific technique of membership determination. In order to 
complete our sample, we have taken the corresponding data 



2 ftp : / /cdsarc . u-strasbg . f r/pub/cats , 
http : / /vizier . u- strasbg . fr 



for these clusters from Loktin et al. (|2004i and the WEBDA 3 
database. The mean proper motions and radial velocities for 
both clusters were computed with data of the ASCC-2.5 and 
CRVAD 4 (Kharchenko et al. I2004at catalogues. All our dis- 
cussions in the following therefore refer to this sample of 652 
clusters. 

Although the current sample contains clusters covering a 
wide range of ages (from logf = 6.6 . . .9.5), it is strongly bi- 
ased towards younger ages. This is not surprising for an un- 
preselected cluster sample: due to dynamical evolution and 
evaporation of members, the typical life time of a cluster is of 
the order of a few hundred Myrs. On the other hand, our cluster 
sample is magnitude limited. Since young clusters contain, in 
general, more luminous stars, they dominate at larger distances 
and their proportion in the sample is somewhat higher than that 
of older clusters. Thus, in average our sample represents the 
moderately young disk. 



3. Spatial distribution in the disk 

Throughout the paper we use the rectangular coordinate sys- 
tem X, Y, Z with origin in the barycentre of the Solar system, 
and axes pointing to the Galactic centre {X), to the direction of 
Galactic rotation (Y), and to the North Galactic pole (Z). The 
corresponding components of the space velocity are U, V, W, 
respectively. Unless other specifications are given, we consider 
the (X, Y) plane to be the Galactic plane. 



3 http://obswww.unige.ch/webda/ 

4 Catalogue of Radial Velocities with 
Astrometric Data can be retrieved from the CDS by 
ftp : //cdsarc . u-strasbg . fr/pub/cats/III/239 



Fig. 2. Open clusters and interstellar extinction. Distributions of the total extinction Ay (left panel), and of the extinction co- 
efficient ay (right panel) in the wider Solar neighbourhood. In both panels the contours show Ay, ay-levels according to the 
grey-colour scales shown at the right upper corners of the panels. The circles represent open clusters. Their sizes (and colour) 
correspond to Ay, ay of the clusters. The blank areas in the periphery of the plots are due to a lack of clusters in our sample at 
these locations. 



3.1. The completeness of the sample and the local 
density enhancement 

In order to estimate the spatial completeness of the sample, we 
studied the distribution of clusters in the (X, Y) plane and com- 
puted their surface density S as a function of the projected dis- 
tance dxy from the Sun. The corresponding histogram is shown 
in the left panel of Fig. [0 Due to the low number of clusters 
at small distances from the Sun, the first bin represents a circle 
with a radius of 150 pc. All other bins are concentric rings of 
100 pc width each. With the exception of a considerable ex- 
cess at dxy = 0.35 . . .0.55 kpc, the distribution is almost flat 
up to = 0.85 kpc. At larger distances, the surface density 
2 of clusters is steadily decreasing. Assuming a uniform den- 
sity model for the distribution of clusters in the Solar neigh- 
bourhood, we can interpret this behaviour as evidence for an 
increasing incompleteness of our sample at projected distances 
larger than 0.85 kpc. We adopt d'L = 0.85 kpc, as the complete- 
ness limit of our sample. 

On the other hand, the enhancement of the observed clus- 
ter density at d„ = 0.35 . . .0.55 kpc is too significant to be 
attributed to random fluctuations in the spatial distribution of 
open clusters in the Solar neighbourhood. Assuming a density 
of 100 + 10 clusters per square kpc as typical of the "clus- 
ter field" around the Sun, we would expect 56 ± 8 clusters at 

= 0.35 . . .0.55 kpc. Instead of this, 91 clusters are counted 
in this region. Does this enhancement describe a general prop- 
erty of the local galactic structure or is this excess caused by a 
particular group (groups) of open clusters with a common evo- 
lutionary history? 



In order to understand the nature of this feature, we con- 
struct the surface density distribution for samples of differ- 
ent ages. Altogether four cluster generations are considered: 
the sample Tl including the 269 youngest clusters with ages 
log? < 7.9; T2 with 101 moderately young clusters of ages 
of log/ = (7.9,8.3]; T3 with 132 intermediate age clusters of 
log? = (8.3,8.6]; T4 with 150 older clusters, logf > 8.6. The 
age ranges are selected in such a way that the pattern of density 
distributions reveals the highest contrast. For each sample we 
construct the histogram of the surface density, smoothed with a 
five-point rectangular filter. The bin size is selected to be 0.025 
kpc. 

The results are shown in the right panels of Fig. [2 The hori- 
zontal dashed lines are average levels of the "field" surface den- 
sity So,/, determined separately for each sample. For the sam- 
ples T2 and T4, they are simple averages of the densities within 
the completeness area. The first bins in both groups, where by 
chance one single nearby cluster produces a maximum of the 
distribution, are not taken into account due to their low statisti- 
cal significance. For the samples Tl and T3 the "field" density 
is estimated from different ranges of distance: Eo,i is computed 
from data in d xy = 0.50 . . .0.85 kpc, whereas for Zo,3 we con- 
sider a distance range of dxy = 0.65 . . .0.85 kpc. We discuss 
these groups in more detail in §|5] 

3.2. Open clusters and interstellar extinction 

Studies of the distribution of the obscuring matter in 
the Galactic disk have been carried out for many years 
(Sharov 1963 FitzGerald [T9681 Lucke Il978l Neckel & 



lable 1. Distribution or the clusters along the Z-axis within a cylinder or radius or d x 
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Klare 119801 Arenou et al. 1992 1. Usually, the methods make 
use of photometric data and spectral classifications of field stars 
over the whole sky and provide both reddening and distance 
evaluation. As a rule, these observations allow a mapping of 
interstellar extinction within the local disk, typically up to 3 
kpc. Recently, Hipparcos parallaxes and data on the equivalent 
widths of Nal lines in spectra of about 1000 local stars were 
used for a detailed mapping of the immediate (d < 0.3 kpc) 
Solar neighbourhood (Lallement et al. 2003 1. As an upper limit 
of extinction the maps constructed in Schlegel et al. ( 1998 ) can 
be taken, although they are of limited use in the 1 kpc neigh- 
bourhood of the Sun. 

Since open clusters reside in the Galactic disk, one must 
expect that they are subject of sometimes strong and irregu- 
lar interstellar extinction which should be taken into account in 
their study. On the other hand, open clusters having a compact 
structure, accurate distances, and with reddening determined, 
as a rule, from a number of stars could be used as suitable 
probes tracing the true distribution of reddening in the Galactic 
disk. The purpose of this section is to reveal the important fea- 
tures of the distribution of interstellar extinction from data on 
reddening of open clusters. 

In Fig.|5]we show the spatial distribution of open clusters 
together with parameters describing the extinction in V: the 
total extinction Ay and the extinction coefficient ay = Ay Id. 
Since Ay characterises the total amount of absorbing matter 
along the line of sight, the total extinction is increasing with 
distance. One observes this increase towards to the periphery of 
Fig.[2](left panel), and notes its different behaviour in different 
directions. Generally, the Ay distribution can be used for ex- 
amining the visibility conditions in a given direction and their 
apparent trends. The extinction coefficient is related to dust and 
can be used for mapping dust clouds (Fig.|2 right panel). 

The total extinction was computed in quadratic areas of 
75x75 pc 2 and then smoothed via a 5x5 rectangular window 
i.e., the pattern displayed is averaged over a scale of about 0.4 
kpc. In contrast to this, for better representation of cloud struc- 
tures, we smoothed the maximum values of ay (instead of their 
averages like for Ay) which were determined in each square of 
75x75 pc 2 . 

In Fig. 13 we show, also, all clusters of our sample with in- 
dividual values of Ay and ay indicated by symbol size (and 
colour). The same size (colour) scale is used in both panels: 
the smallest (white) circles mark open clusters with Ay < 0.5 
mag (left panel) or ay < 0.5 mag/kpc (right panel); the small 
(yellow), medium (dark yellow), large (red) and the largest 
(brown, dotted) circles indicate clusters with Ay/ay ranging 



into (0.5, 1.5], (1.5,2.5], (2.5,3.5], and > 3.5 mag or mag/kpc, 
respectively. 

Some general features of the extinction pattern can be ob- 
served in the left panel of Fig. 13 At the location of the Sun we 
find a cavity with a small extinction of the order of 0^1 . . . 0^2, 
which is transformed into a long transparent (Ay m 0? 1 2 . . . 0?3) 
corridor located in the directions I m 80° and 220° of galac- 
tic longitude. Another corridor of about the same transparency 
(Ay * 07"2) could be seen if we consider only clusters with 
Z < -20 pc. It is directed towards / * 120° and I * 300° from 
the Sun. 

Another feature is the complete absence of clusters (and 
also extinction data) in some more peripheral areas. Most 
interesting are the areas at (X, Y) w (-1.2,1) kpc and at 
(X, Y) « (1.2, l)kpc which manifests a screening effect of mas- 
sive nearby clouds. Some of these clouds are clearly seen in 
the right panel, and they are identified with the Auriga-Taurus 
clouds at X < -0.1 kpc, and the Ophiuchus clouds at X > 
kpc. However, the empty area at (X, Y) m (1.2, 1) kpc is not 
caused by the Ophiuchus clouds. Instead, the molecular clouds 
located in the Aquila Rift region play a more important role 
here. The clouds are extended farther in the direction of I w 50° 
(Dame et al. 2001), and are not seen in our map because the 
known clusters in this direction are too faint in the optical to be 
identified in the ASCC-2.5. 

The impact of the clouds on the visible distribution of the 
clusters is not limited to the Solar vicinity, only. Their shadow 
can be seen at larger distances and they influence locally the 
visibility of such structures as the Sagittarius and Perseus arms 
at some kpc away from the Sun (see § 13. 4> . Although the spatial 
density of open clusters is lower than the density of individual 
stars, the extinction pattern obtained is in a good agreement 
with most of the previous publications based on data of field 
stars. 

3.3. Distribution of clusters perpendicular to the 
Galactic plane 

Distribution along the Z-axis is usually described in terms of 
the barometric formula given in eq. where D(Zq) is the den- 
sity at the position of the symmetry plane Zo and hz is the scale 
height: 

D(Z) = D(Z )exp|-^^J. (1) 

The distribution parameters Zo, hz and Z)(Zq) were computed 
with data on 259 clusters located within the completeness limit 



Fig. 3. Distribution of clusters perpendicular to the Galactic 
plane. The distributions are constructed for clusters within the 
completeness limit of 850 pc. Panel A: all clusters, panels Tl 
- T4: the 4 age samples (from Tl to T4) defined in § 13.11 The 
thin black curves indicate smoothed spatial densities. The thin 
(magenta) curves in panels (T1),(T3) and (A) mark "field" sub- 
samples of the corresponding ages. The vertical solid lines give 
positions (Zo) of symmetry planes, and the dashed vertical lines 
show Zo of "field" subsamples. The thick solid curves cor- 
respond to Gaussian distributions computed with parameters 
listed in Tabled Thick long-dashed lines give the barometric 
distributions fitted to observations. 

determined in § 13.11 Additionally, we considered the distribu- 
tions of clusters of the age groups defined in § 13.11 Further, we 
separated "field" subsamples i.e., excluded probable kinematic 
members of the cluster complexes discussed in §[5] The corre- 
sponding distributions along the Z-axis were constructed with a 
step of 1 pc and then smoothed by use of a rectangular 15 -point 
filter. The results are shown in Fig [3] 

Due to poor statistics (the number of clusters in the different 
groups varies from 45 to 259) and large density fluctuations, a 
simultaneous solution of eq. did not yield a stable and ev- 
ident result. Therefore, we computed the distribution parame- 
ters in two steps. At first, the position of the symmetry plane 
was assumed to be a simple average in each group. The me- 
dians and modes were also considered but the results did not 
differ significantly. For each sample, the standard deviation <x z 
of |Z - ZqI was determined as a statistical measure of cluster 
scattering. In a second step, the parameters D{Zq) and hz were 
derived by use of the integrated form of eq. 0: 



where is the surface density of clusters in a layer with 
width 2 |Z - Zo|. The parameters were found from a non-linear 
fit of eq. to the observed distributions. The results with cor- 
responding rms errors are given in Tabled 

According to Fig. [3] and Table ^ open clusters form an 
extremely flat system in the Solar neighbourhood. The distri- 
bution demonstrates non-regular behaviour with many peaks 
and gaps. We attribute this to a possible impact of interstel- 
lar clouds on the apparent cluster distribution (see also § I3.2> . 
Although varying from one sample to another, the symmetry 
plane of the cluster system is located below the {X, Y) plane. 
This finding indicates the position of the Sun above the physi- 
cal plane of the system of open clusters. The derived displace- 
ment Zo = 22 + 4 pc is in good agreement with recent re- 
sults obtained from studies of other young Population I ob- 
jects: according to Reed (2005) the position of the symmetry 
plane of OB stars within 1200 pc from the Sun is located at 
Zo = 19.5 ± 2.2 pc, whereas the symmetry plane of the dust 
layer is at Zo = 22.8 ± 3.3 pc (Joshi 2005 1. The scale heights 
vary within 50-60 pc without an evident trend with respect 
to the ages of clusters (the youngest (Tl) and the oldest (T4) 
groups have the largest hz value). Although the standard devi- 
ations cr z are always larger than the scale heights, the corre- 
sponding Gaussian and barometric distributions do not differ 
significantly (see Fig.[3J. 

We found the total volume and surface densities of clusters 
in the plane of symmetry to be D{Zq) = 1015 kpc~ 3 and S = 
1 14 kpc 2 , respectively. The latter number considerably (by a 
factor of about 5) exceeds the surface density of open clusters 
which is inferred as £ =24.9 kpc~ 2 from Janes et al. (1988 
Table VIII). This discrepancy can be explained by mainly two 
reasons: a) a significant number of clusters within 0.85 kpc in 
our sample is missing in the sample of Janes et al. {I988l, and 
b) the completeness of their cluster sample is overestimated 
by Janes et al. (I1988> . A combination of these effects should 
increase their result for the surface density by a factor of 3-6. 

Provided that the cluster density distribution in the Galactic 
disk can be described by a radial density profile with a scale 
length of the disk hi, we estimate the total number of clusters 
N tot in the Galaxy as 




Here Rq is the distance from the centre of the Galaxy, R ™ 1 
is the radius of the Galactic disk, and hi is the scale length, 
adopted from Bahcall & Soneira ( 1980 > to be hi - 3.5 kpc. 
With Rq =15 kpc, we derive the total number of open clus- 
ters in the Galaxy as N, nt = 93 000. Within a reasonable range 
of the parameters R l ™ = 13-20 kpc, and hi > 3 kpc the 
above estimate is valid within 20%. For smaller values of hi, 
the expected number of open clusters would come out to be 
considerably larger. 



Fig. 4. Distribution of open clusters in the .XT-plane centered at the location of the Sun. Blue squares mark Tl clusters, black 
plusses - T2, red circles - T3, and black crosses - T4 clusters. Dots indicate the positions of clusters with known distances and 
ages from the DLAM catalogue (v2.0), which are not included in our sample. The youngest (log t < 7.9) of these clusters are 
marked by triangles. Long dashed and solid thick curves are the grand design spiral arms with a pitch-angle of -6° fitted to 
young clusters in the Perseus and Carina-Sagittarius regions. The large circle around the Sun indicates the completeness area of 
our sample. 



3.4. Distribution of open clusters in the Galactic plane 

The distribution of open clusters in the Galactic plane is pre- 
sented in Fig.|4] There we show also some 300 clusters from the 
catalog DLAM (version 2.0) which have distances and ages, 
but are missing in our sample since they are mostly too faint to 
be identified in the ASCC-2.5 data. Only about 30 of them are 
located within the completeness area of our sample. It turned 
out that all these 30 clusters did not fulfill our criteria for being 
real open clusters and were thus excluded from the sample. The 
retained clusters shown in Fig.@]append the current distribution 
at larger distances. 

Until now the pattern of the spiral structure of the Galaxy 
is not clearly established. In order to fit the observations, dif- 
ferent models have been proposed which suggest two, three or 
four spiral arms, with pitch angles ranging from -5° to -21° 
(see Valee 1995 for reference). Also, many attempts were un- 
dertaken to discuss the apparent distribution of open clusters in 
the Galactic plane with their relation to a local spiral pattern 



(see e.g. Becker [T963l Janes et al. H988i . Since the clusters of 
our sample have homogeneous distances and ages, we can try 
to fit their distribution by a spiral pattern, but we must take into 
account that the real distribution of open clusters in the optical 
is hidden by lumpy extinction, and we see distant clusters only 
by chance through transparent windows (cf. Figs.|2]and|4j. 

In Fig. |4] we show the grand design pattern of logarith- 
mic spirals with a pitch angle p = -6° fitted to the Perseus 
and Carina-Sagittarius arms. The choice of the pitch-angle is 
dictated by our intention to provide a reasonable fit of a sim- 
ple pattern (i.e., a two-arm spiral) to the observed distribu- 
tion of young clusters at different galactocentric distances. We 
have succeeded in getting a reasonable agreement along the 
whole Perseus arm (from I w 90° to about 270°), the Carina- 
Sagittarius arm (from I w 300° to about 360°), and its extension 
in Vulpecula at about (X, Y) — (1, 1.5) kpc. Taking into account 
relative errors of distances of about 15%, it is not excluded that 
a few distant clusters at about X=3 kpc and X=5 kpc trace the 



lable 2. Kinematic and orbital parameters or the sub-system or 
open clusters in the Galaxy 



Parameter Value rms error 



Uq, km/s 


+9.44 


1.14 


V G , km/s 


+ 11.90 


0.72 


W e , km/s 


+7.20 


0.42 


o~u, km/s 


13.86 


0.81 


o>, km/s 


8.75 


0.51 


o"h/, km/s 


5.05 


0.30 


cry : cr v : cr w 


1 : 0.63 : 0.36 




R a , kpc 


8.631 


0.034 


c* > kpc 


0.427 


0.024 


R p , kpc 


6.706 


0.056 


kpc 


0.682 


0.040 


^rntai kpc 


0.260 


0.016 


^ kpc 


0.189 


0.011 


e 


0.127 


0.003 


cr e 


0.037 


0.002 


A, km/s/kpc 


+ 14.5 


0.8 


fi, km/s/kpc 


-13.0 


1.1 



inner convolutions of the Perseus and Carina-Sagittarius spi- 
rals. 

No open clusters are visible at the expected position of the 
Carina-Sagittarius arm at / = 20 - 55°. We explain this by mas- 
sive molecular clouds located in this direction (see the discus- 
sion before). According to Bica et al. (2003b I the gap is well 
filled by 45 embedded clusters detected in the 2MASS atlas. 
However, they cannot be used as tracers of the spiral arm since 
reliable distances are still unknown for them. 

4. Kinematics of open clusters 

The first efforts to transform known proper motions of open 
clusters to a common reference system was undertaken by van 
Schewick (I197U who derived absolute proper motions of about 
70 open clusters in the FK4 system. Lynga ( 1982 1 provided ra- 
dial velocities for 108 open clusters, but only for a few tens of 
them all three vectors of the space velocity could be obtained 
at that time. Due to the Hipparcos observations and to the de- 
velopment of coravel-type spectrovelocimeters, this situation 
has improved considerably. Currently, all 652 clusters of our 
sample have accurate proper motions in the Hipparcos system 
and uniform heliocentric distances. For about 55% of them, the 
mean radial velocities (RV) are also available (Paper II). Since 
the kinematic data of clusters are generally more accurate than 
these of single field stars, we make use of this advantage to 
learn about the kinematics of the Galactic disk, in more detail. 
On the other hand, we should keep in mind that the subsample 
of open clusters with complete kinematic data is biased towards 
younger clusters: RVs are available for about 75% of clusters 
younger than logf = 8.3, whereas for older clusters, radial ve- 
locities are known for only 3 1 % of them. This has to be taken 
into account if one compares the kinematic parameters derived 
from proper motions (or tangential velocities) of 652 clusters or 
from the subsample of clusters with proper motions and radial 
velocities. 



In the following, we distinguish three types of space veloc- 
ities. Firstly, these are observed velocities which are computed 
directly from COCD proper motions and radial velocities. For 
each cluster, this space velocity contains systematic compo- 
nents due to Galactic differential rotation and due to Solar mo- 
tion with respect to the centroid of the open cluster system, 
as well as a peculiar motion of a cluster with respect to this 
centroid. The second type are velocities corrected for Galactic 
differential rotation, and, finally, the third type are velocities 
corrected both for Galactic differential rotation and for Solar 
motion (i.e., the peculiar velocities of clusters with respect to 
the centroid of the open cluster system). 

As a first step, we determine kinematic parameters describ- 
ing the basic motions of the system of open clusters in the 
Galactic disk: the motion with respect to the Sun, and the dif- 
ferential rotation around the Galactic centre. The Solar motion 
components U ,V e , W© were derived with the complete kine- 
matic data of 148 clusters located within d'L < 0.85 pc as: 

U Q = -U, V = -V, W = -W, (4) 

with U, V, W computed from proper motions and radial veloci- 
ties with 

U = k d (fii sin / + nt cos / sin b) — RV cos / cos b, 
V = kd(-fii cos/ + fib sin/ sin/?) - RV sin/ cos/?, (5) 
W = -kdfib cosb-RVsinb, 

where k is the scale factor converting angular motions into lin- 
ear velocities. The Oort's constants A and B of the Galactic 
rotation were determined from the proper motions /ji, fib of 58 1 
clusters with distances less than d xy < 2500 pc by a LSQ- 
solution of the equations: 

kdfii = (U sin/ - V Q cos I) + (A cos 2/ + B) cos b, 
kdfib = (Uq sinfo cos/ + V sin/? sin/- W cosb) - 

— sin 2b sin 2/. (6) 

A rigorous solution of eq. (|5j assumes velocity components 
corrected for Galactic rotation, whereas the velocities in eq. 
should be free from the Solar motion. In order to fit these in- 
troduced corrections separately, we solved the equation system 
by iteration. Only a few iterations were necessary before 
convergence was achieved. 

In the next step, the space velocities of 148 clusters within 
d xy < 0.85 kpc corrected for Solar motion and Galactic rotation 
were used for computing the parameters {cr u , cr v , <j w ) of the 
velocity ellipsoid. Assuming the form of the Galactic potential 
proposed by Saio & Yoshii dl979> . we also derived elements 
of their Galactic box orbits. This includes the apocentre and 
pericentre distances R a and R p , the eccentricity e and the max- 
imum vertical distance Z max which a cluster can reach in its 
orbital motion. The results are given in Tabled 

Fig-0shows the distribution of kinematic data of open clus- 
ters versus Galactic longitude for different ranges of distances 
as indicated at the top of the figure. The curves mark system- 
atic contributions of Solar motion and Galactic rotation to the 



Fig. 5. Observed kinematic data of open clusters versus Galactic longitude for four ranges of distances (as indicated at the 
top). The upper panels are for proper motions, the middle - for tangential velocities, and the bottom - for radial velocities. The 
parameters of individual clusters are marked by thin horizontal strokes together with the vertical error bars (in the RV panel the 
error bars are given only for clusters for which the rats-errors could be estimated). The crosses in the rightmost panels indicate 
clusters with d > 3.5 kpc. For a given distance range, the curves show systematic velocity components due to Solar motion 
and Galactic rotation computed with the parameters £/©, V , W Q , A and B in Table|2] The thick dashed curves correspond to the 
contributions at the smallest distance of a given distance range, whereas the thin solid curves show this contribution at the largest 
distance. For example, in the upper left panel the dashed curve shows systematic variations versus longitude in proper motions 
expected for a cluster at 0.085 kpc from the Sun, whereas the solid curve shows this variations for a cluster at 0.5 kpc. 



observed velocities of open clusters which are derived with the 
parameters in Table 13 Since this systematic velocity compo- 
nent is a function of distances, we show its contribution with 
the dashed and solid curves for the minimum and maximum 
distances in a given distance range, respectively. The rms er- 
rors for proper motions and radial velocities are taken from the 
COCD, the errors of tangential velocities are derived as: 

4=tfd 2 (i 2 (6l + 6% (7) 

where 5^ = s^f/u, 5d = sj/d are the relative errors of a proper 
motion and distance. The relative error of a proper motion 5^ 
is computed from the individual data of clusters, whereas Sd is 
assumed to be 15% of a distance (see Paper II). 

In an ideal case, the observed velocities should be located 
between the dashed and solid curves in each panel of Fig. [5] 
Deviations significantly larger than the individual rms errors 
indicate peculiarities in cluster motions. For example, a few 
clusters at longitudes 44. . . 1 24° and distances 0.5 ... 1 .5 kpc show 
unusually large proper motions and tangential velocities. We 
found that all of them are older than log t = 8.6. With the 



exception of NGC 6991 their radial velocities are unknown, 
so we cannot determine the complete space velocity vector. 
Fortunately, their location close to / = 90° allows us to esti- 
mate the U -components, and we found that these clusters move 
towards the Galactic anticentre: U -velocities corrected for the 
Galactic rotation are from -30 to -50 km/s. Direct calculations 
of the velocity vector for NGC 6991 give a similar radial com- 
ponent U = -58 km/s. The age of the clusters as well as their 
motion towards the Galactic anticentre remind us of the Hyades 
(log t = 8.90, U = -43.5 km/s). We consider this cluster group 
in §|5]in more detail. 

Oort's constants of Galactic rotation are related to the dis- 
tance of a cluster from the Galactic centre Re, and to the angu- 
lar velocity Q. of the disk as: 

* - -mi 

B = A - Q 

Assuming the galactocentric distance of the Sun /?c,o = 8.5 kpc, 
we obtain the angular velocity of the Galactic rotation Qq = 



A - B — 27.5 + 1.3 km/s/kpc and the corresponding rotation 
velocity of the system of open clusters as 233.8 km/s. 

This result is in a good agreement with the rotation param- 
eters derived from Hipparcos data over a similar range of ages 
(less than 1 Gys). For example, Feast & Whitelock ( 1997 1 ob- 
tained (A, B) = (14.82 + 0.84,-12.37 ± 0.64) km/s/kpc from 
kinematic data on 220 Galactic cepheids located mainly within 
4 kpc from the Sun. Using data on about 2000 O-B stars, 
Torra et al. (2000 1 found (A,B) = (11.8 ± 0.4,-12.3 + 0.4) 
km/s/kpc for d < 2 kpc. From data on about 1820 0-B5 stars 
within 1 kpc from the Sun, Branham (2002 1 derived (A,B) = 
(+14.9 ± 0.8,-15.5 ± 0.7) km/s/kpc. 

Also, the parameters of Solar motion derived with open 
clusters agree quite well with previous results for the moder- 
ately young disk. Based on a kinematically unbiased sample 
from the Hipparcos catalogue in the Solar velocity, Dehnen & 
Binney ( 1998 1 determined the Solar motion with respect to the 
LSR as(t/,V,W) = (10.00 ± 0.36,5.25 + 0.62,7.17 + 0.38) 
km/s. Whereas the U and W-components are in a perfect agree- 
ment with the results derived with open clusters, the rotation 
component is about two times smaller than in Table |2] The 
difference is mainly caused by the different evolution param- 
eters of the samples used. The sample of Dehnen & Binney 
contains Main Sequence stars from the Solar neighbourhood 
with (B — V) < 1.2 i.e., stars of different ages, from a few 
million years up to 10 Gyr. As we already mentioned before, 
our cluster sample represents the population of the moder- 
ately young disk, (see also § |6j. Indeed, comparing our re- 
sults with that of Torra et al. (2000 1 who found (U, V, W) Q = 
(11.00 ± 0.2, 12.9 + 0.2,6.8 + 0.1) km/s with OB stars, we see 
a good agreement for all components of the Solar motion. We 
conclude that both systems of young stars and of open clusters 
move similarly. 

Since the age for each open cluster of our sample was deter- 
mined by the same method, we can check whether real trends 
in kinematic parameters exist, depending on ages. The most in- 
teresting trend is a dependence of the cosmic (or true) velocity 
dispersion on the cluster age. Due to the inhomogeneity of the 
cluster RVs, it is rather difficult to estimate a realistic contri- 
bution of the rms errors in the measurements to the observed 
dispersion of space velocities. Therefore, we prefer to use the 
dispersion of tangential velocities which can be properly cor- 
rected for the rms error of observations. Let ej be the mean 
error of tangential velocities of a given cluster sample, then 

0~j — CTj- — Sj (9) 

where &j is the wanted cosmic dispersion, and cr T is the dis- 
persion of observed tangential velocities corrected for the Solar 
motion and Galactic rotation. We consider all clusters within 
d c xy , and with logf from 6.7 to 9.5 binned in 10 overlapping 
groups of ages. The resulting relation is well approximated by 
the following equations: 

d- Tl = (1.44 + 0.45) log f 6 + (3.93 + 0.84), (10) 
& Th = (1.33 ±0.17) log t 6 + (1.65 ±0.33) 

where is the cluster age in units of Myr. For a given log f, 
(<5"r,, <5y 6 ) represent the expected dispersions of tangential ve- 
locities of clusters of this age and describe the ellipse of the 



velocity dispersion, or simply the velocity ellipse. Since the 1 / 
components might be biased by residual rotations, we use the 
second equation of dlOt for estimating the ellipsoid of the cos- 
mic dispersion. Moreover, for a flat system of open clusters the 
Ti, component does not differ significantly from W. Using the 
ratio (cry : cry '■ cr w ) in Table |5J we derive the ellipsoid of 
the cosmic dispersion as (7.2 : 4.1 : 2.6) km/s for young open 
clusters with logf = 6.7, and (17.5 : 10.0 : 6.3) km/s for old 
clusters with log t — 9.5. 

The derived relations were compared with recent findings. 
Dehnen & Binney ( 1998 1 determined dispersions for MS stars 
in different colour ranges. Their youngest group was selected in 
the (B - V)-range from -0.238 to 0. 139. According to our esti- 
mates, the average log t of these stars should be somewhere be- 
tween 8.2 and 8.8. Then, fromeq. dlOt we derive d~w = 4.6. ..5.4 
km/s. The corresponding value of Dehnen & Binney ( 1998 1 is 
5.46 ± 0.61 km/s and fits well to our estimates. According to 
Torra et al. ( 2000 1, the velocity dispersion of OB stars increases 
from 4.3 km/s for about 15 Myr old stars to 5.5 km/s for stars 
older than 120 Myr. These results are comparable to ours, al- 
though they are larger by about 20%, than those derived from 
eq. E5- 



5. The local density inhomogeneities 

In this section we discuss the density enhancement of open 
clusters found in the spatial distribution in the Solar neighbour- 
hood illustrated in Fig.^ The density peaks are detected nei- 
ther only in one age group (this would be the case if a peak is a 
signature of one physical group of clusters), nor in all of them 
(which would indicate a general property of the cluster distri- 
bution). On the contrary, the clusters from groups T2 and T4 
show distributions which are compatible with a constant den- 
sity of clusters, whereas the samples Tl and T3 have signifi- 
cant maxima within the completeness area. Moreover, the sam- 
ples Tl and T3 display differences in their density distribution. 
In addition to the prominent peak at about 0.4 kpc, the sam- 
ple Tl shows a secondary density excess of lesser significance 
at lower distances which can be interpreted as evidence for a 
prolate structure. The sample T3 demonstrates only one, but a 
broader density excess. In the following we consider the prop- 
erties of each sample separately, and try to answer the question: 
are these density peaks manifestations of independent "cluster- 
ings" of open clusters? If so, to which degree are they separated 
from each other, and from the general cluster substratum? 

Hereafter, we call such cluster groups "open cluster com- 
plexes" (OCCs), in analogy to large star forming complexes 
which include associations and young clusters (see e.g., Janes 
et al l 19881 Efremov fT995l . We would like to stress that unlike 
star forming complexes, which are considered as aggregates of 
objects related to star formation, the OCCs have been primarily 
identified within a sample of classical open clusters. Of course, 
some of them could (but need not) coincide with star form- 
ing regions. In the following, we refer to clusters, which are 
not found to be members of OCCs, as "field" clusters of the 
Galactic disk. 



Fig. 6. Vector point diagrams of tangential velocities of open 
clusters in four age samples (from Tl to T4). The velocities 
are corrected for Solar motion and differential rotation of the 
Galaxy. Dots are clusters (independent of age) within the com- 
pleteness area. The crosses mark clusters of a given age sam- 
ple. The circles and squares are clusters located within the areas 
of the density excesses (samples Tl and T3). Their sizes indi- 
cate kinematic membership probability to belong to an OCC 
correspondingly, the largest symbols indicate lo - members, the 
intermediate circles are 2<x members, and the smallest circles 
are 3<x members. The squares (in T3) mark the members of 
the Perseus-Auriga group (see § I5.3> . The triangles (in T4) in- 
dicate possible members of the Hyades cluster group (see §0] 
and § I5.4> . The plusses in the bottom-left corners in each panel 
give the mean errors of tangential velocities of open clusters. 
The (green) cross-hair indicates the average tangential motion 
of each group, whereas the ellipses show the velocity ellipses 
of the complexes. 



5.1. Tangential velocity and membership in complexes 

In order to get additional pros and/or contras for the physical 
nature of OCCs, we consider the distribution of kinematic pa- 
rameters of open clusters. Since radial velocities have not been 
measured for all clusters, we restrict ourselves to the study of 
two-dimensional kinematics only. As the clusters are located at 
different distances from the Sun, we use tangential velocities 
instead of proper motions. 

In Fig.[6]we compare the distributions of clusters of differ- 
ent ages in the (Tt,, Tf) diagrams (i.e., the vector point diagram 
of tangential velocities, the VPD). The T\,, Ti velocities are cor- 
rected for Solar motion and differential rotation of the Galaxy. 
Here we consider only clusters located within the completeness 
area (dxy < 0.85 kpc) and compare their distribution with dis- 
tribution of clusters in the enhancement areas (in the samples 
Tl and T3). 

A proper selection of kinematic members requires objec- 
tive criteria taking into account peculiarities of the kinematic 
data used. Although we may apply a technique of probabil- 
ity determination normally used for the selection of members 



of a cluster, we must introduce some principal modifications. 
Determining the members from the proper motion distribution 
within a small cluster area, we may assume that proper motions 
of cluster members differ from the average only by their rms 
errors, and generally, the internal motion is negligibly small 
compared to uncertainties of proper motions. This assumption 
is not always true if we consider the common motion of clus- 
ters belonging to an OCC. According to eq. i ll Oi l, we would 
expect a cosmic dispersion comparable to the rms errors of the 
tangential velocities. Due to the definition of the tangential ve- 
locities, even clusters with similar space velocities should show 
different tangential velocities as a function of their galactic co- 
ordinates. This effect will introduce an additional spread of data 
points in the VPD. Since the clusters of our sample represent a 
flat system, this spread should be more significant in T\ than in 
Tb- Let us assume an ideal (i.e., without internal rotation and 
velocity dispersion) stream of clusters located in the Galactic 
plane and well distributed around the Sun. Then, the depen- 
dence of Ti on galactic longitude I for each cluster has a sine 
function: 



T, = -A Tl sin(Z - l ) 



(11) 



where lo is the apex direction of the stream and At, is the max- 
imum contribution of the velocity component T\ to the space 
velocity of the stream. An additional rotation will simply make 
the mean tangential velocity T\ of a cluster group different from 
zero. It is self-evident that the spread of tangential velocities in 
the VPD will be smaller if OCC members are located at similar 
galactic longitudes as well as when they have low common ve- 
locities with respect to the centroid of open clusters. However, 
in reality the dependence will be more or less disturbed by a 
spread in the location of clusters along the Z-axis, by possible 
internal motions, uncertainties of kinematic data, and so on. 
Below we discuss the kinematic behaviour of the complexes in 
further detail. 

Due to the reasons described above, we cannot consider the 
differences of cluster motions with respect to their group aver- 
age to be a measure of the membership probability as we did in 
the case of stars in a cluster. In the case of OCCs, we must con- 
sider residual motions of individual clusters with respect to the 
velocity ellipse defined by a realistic dispersion of tangential 
velocities. 

The membership probability P, of a cluster to belong to a 
OCC is computed as 



P t = exp 



V Ti I 



\ Tl, J 



(12) 



Here s' T are rms errors of tangential velocities of the z'th clus- 
ter, and ST' are residual motions of the cluster with respect to 
the velocity ellipse of the OCC considered. Depending on the 
residual motions 5T', we classify the membership probability 
as high (dr/s'j. < 1, P, > 61%), moderate (1 < ST/dj, < 
2, P, = 14% - 61%), and low (2 < SP/d T < 3, P, > 1% - 
14%). Traditionally, we call such objects lcr, 2<x, and 3<x mem- 
bers (see Paper I), respectively. Note that all clusters within 
the velocity ellipse have the highest membership probability 
(P, = 100%). 



1 he parameters or velocity ellipses were chosen in accor- 
dance with equations dlOt derived in § [4] The centres of the 
distribution and membership probabilities were computed in an 
iterative procedure. As a first approximation of the centres, we 
used a simple average of the tangential velocities of all clus- 
ters at a given age and we derived membership probabilities 
by eq. H2i . After that, we selected the lcr-members to derive 
the next iteration of the distribution centre. The iterations were 
stopped when the list of the lcr-members did not change in sub- 
sequent iterations. The results are shown in Fig.[6]and Table[7] 

5.2. OCC 1: a complex of clusters associated with 
Gould's Belt 

For the sample Tl, we found 37 clusters within d xy = 0.5 kpc, 
with 28 of them located within the main maximum of the den- 
sity distribution at d_ n , = 0.325 . . .0.5 kpc. Recomputing the 
surface density distribution in terms of numbers of clusters, 
we infer that 16±4 clusters should belong to the field, whereas 
20+4 clusters are responsible for the density excess. Thus, we 
may expect about 20 members in a potential OCC. About the 
same number of clusters is also obtained from the statistics of 
kinematic membership. Table[7]shows that within 0.5 kpc from 
the Sun, we found N = 23 clusters which are lcr members from 
the distribution of their tangential velocities. We conclude that 
the lcr threshold for kinematic probability (P, =61%) sepa- 
rates well the complex members from the field. 

The distribution of the OCC 1 candidates in 3D-space is 
shown in Fig. [7] For convenience, we present the results in a 
coordinate system related to this cluster group. The coordinate 
system X', Y', Z is rotated in the Galactic plane to such that the 
X'-axis corresponds to the line of nodes of the OCC 1 plane. 
The centre of the system is shifted to the apparent centre of the 
group by AX = -0.078 kpc and AY = -0.053 kpc. 

According to Fig. [7] the OCC 1 clusters constitute an elon- 
gated configuration immersing below the galactic plane in the 
region of the Orion association by about 200 pc, and then rising 
above the galactic plane by about 100 pc at the opposite direc- 
tion in the IC 4665 region. The apparent centre of the com- 
plex is located close to the Sco-Cen group with the clusters 
Mamajek 1 and Platais 8. On average, the OCC 1 resides be- 
low the Galactic plane (Z = -46 pc for the lcr members). The 
ascending node of the OCC 1 plane is located near the cluster 
Alessi 5 at (l,b) = (288°08, -2°05), whereas the descending 
node is somewhere near ASCC 127 at (l,b) = (112°26,4°13). 
This finding is well compatible with the results by Torra et 
al. (|2000 1 who derived the longitude of the ascending node of 
Gould's Belt as I = 275 . . . 295° from the Hipparcos data of 
OB field stars. It seems that the density excess of Tl clusters 
can be associated with the Gould's Belt complex which is usu- 
ally believed to contain early-type field stars, associations, or 
interstellar gas (e.g., see a review by Poppei n"997> . 

The OCC 1 has a flat structure, but is not strictly planar. 
In the (Z, T')-pl ane it shows a significant inclination (as the re- 
gression solution shows) of 19° with respect to the Galactic 
plane. Torra et al. (2000) found an inclination of ; = 16 . . . 22° 
from the data on field stars. In the (Z, X')-plane, the complex 



lable 3. Candidates of the OCC 1 (the Gould s Belt complex) 



COCD 


Name 


X 


Y 


Z 


log? 


P, 






kpc 


kpc 


kpc 


yrs 


% 


41 


Stock 23 


-0.291 


0.244 


0.014 


7.51 


25 


42 


Melotte 20 


-0.159 


0.102 


-0.021 


7.55 


2 


46 


IC 348 


-0.354 


0.125 


-0.120 


7.79 


84 


68 


Collinder 65 


-0.301 


-0.047 


-0.058 


7.41 


100 


72 


Collinder 69 


-0.413 


-0.113 


-0.092 


6.76 


100 


73 


NGC 1981 


-0.334 


-0.178 


-0.130 


7.50 


100 


74 


NGC 1976 


-0.329 


-0.183 


-0.132 


7.71 


100 


75 


NGC 1977 


-0.415 
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presents a less pronounced structure. It is a rather smooth tran- 
sition from the lowest location near NGC 1980 to Stephenson 1 
via Vela OBI and Mamajek 1, and then a steeper rising from 
IC 348 to IC 4665 via ASCC 127. 

In Table we list all Tl clusters located within d xy = 0.5 
kpc from the Sun. Due to its high membership probability (P, = 
100%), we extended this list by one Orion cluster (NGC 1980) 
residing at d xy = 5 1 8 pc. For each cluster in Table|3] we provide 
the COCD number, the cluster designation, the (X, Y, Z) coor- 
dinates, the cluster age, and the kinematic membership proba- 
bility P, for belonging to the Gould's Belt complex. 

As described above, clusters located within the velocity el- 
lipse got the highest membership probability P, = 100%. On 
the other hand, some clusters like e.g. Melotte 20 (a Per) which 
were believed to be members of the Gould's Belt system (e.g., 
see Olano 1982 1 have a low kinematic probability. 



Fig. 7. Spatial distribution of the OCC 1 clusters in 3D-space. The proper coordinate system X', Y',Z of the OCC 1 is described 
in the text. The coordinate axes X and Y of the Galactic rectangular system are shown with dotted lines. Two circles centered 
at the origin of the Galactic coordinate system indicate the outer boundary of the area of the enhanced density at 0.5 kpc and 
the completeness area at 0.85 kpc. The solid ellipse outlines the OCC 1 members selected as described in the text. The dashed 
ellipse represents the Olano ( 1982 1 model of gas in Gould's Belt, given here for comparison. The small dots are all the clusters 
from our sample, the crosses mark Tl clusters, and the filled circles are for member candidates of the OCC 1. The small (cyan) 
circles are 3cr members, the dotted circles (green) of medium size show 2<x members, and the largest dotted circles (light yellow) 
with a (green) core indicate lcr members. The grey colour contours give the projected density distribution of Tl clusters, derived 
separately in each panel. The density distribution was computed in windows of 50 x 50 pc 2 and smoothed with a 5-point rect- 
angular filter. The rules to the right indicate the corresponding density scales where the (X',Z) and (Z, Y') planes share the same 
grey colour scale. The straight line in the panel (Z, Y') is the regression line computed with lcr clusters. The resulting inclination 
i is given in the left upper corner of (Z, Y') plane. 



Fortunately, all the lcr members of OCC 1 have full space 
velocity vectors. Therefore, we can analyse the spatial mo- 
tion of this complex in more detail. On average, OCC 1 ro- 
tates around the Galactic centre at about the same velocity as 
the centroid of the open cluster system (the relative velocity 



is V — -0.2 +1.1 km/s). The average tangential velocity of 
the most probable members of the OCC 1 is Ti = -3.6 ±1.4 
km/s i.e., it differs significantly from zero. Therefore, we may 
conclude an internal rotation of the OCC 1 complex around its 
centre near IC 2391 and Mamajek 1, and in the direction of 



Galactic rotation, bxcluding these two clusters which do not 
contribute to the determination of the internal rotation, we de- 
rive T/ = -2.9 ± 0.7 km/s. At an average distance of 0.39 kpc 
from the OCC 1 centre, one obtains the corresponding angular 
velocity as co — -8.6±2.2 km/s/kpc. This is considerably lower 
(in absolute units) than the value a> m -20 km/s/kpc derived by 
Bobylev (2004) from the kinematics of Hipparcos stars pro- 
posed to be members of associations and open clusters. 

The expansion velocity of the OCC 1 computed as a resid- 
ual motion along a radius-vector projected onto the XY plane 
and with origin in the OCC 1 centre is V r = 2.1 + 1.1 km/s. This 
seems to be insignificant. The latter conclusion confirms the re- 
sults of Torra et al. (2000 1 who found an expansion velocity of 
1.5 ± 1 km/s for young field stars at d > 0.25 kpc. 

The total velocity dispersion cruvw for the 23 members of 
the OCC 1 

v 2 uvw = (°u + cr v + c^)/ 3 

is found to be 4.8 ± 0.7 km/s. This is significantly less than 
the total dispersion of 7.6 + 0.7 km/s computed from all 69 
Tl clusters in the completeness area. This difference gives an 
additional confirmation for the distinction of the complex from 
the general field. 

5.3. OCC 2 and the Perseus-Auriga group: two in one 

According to the cluster counts within d„ = 0.65 kpc, the sam- 
ple T3 contains 53 clusters, with 48 of them in the location of 
the main density maximum (i.e., at d„ = 0.325. ..0.65 kpc). 
The lower density of T3 clusters at d xy < 0.325 kpc from the 
Sun can well be seen in Figs.[0and|8] Recomputing the surface 
density distribution in terms of numbers of clusters, we infer 
that 33±6 of them belong to field clusters and 20±6 form the 
density excess. The number of kinematic candidates of OCC 2 
(i.e., lcr members) is 27 (see Table0. Although a difference of 
7 clusters is still within the statistical uncertainty, the disagree- 
ment between both estimates is larger than in the OCC 1 case. 
The average density used for computing the number of field 
clusters is determined from the distribution of T3 clusters up to 
= 0.85 kpc, and it is slightly higher than the density in the 
inner part of d xy < 0.325 kpc from the Sun (Fig.Q. Assuming 
the density in the Solar neighbourhood as being representative 
for the T3 group, we obtain 20 field clusters and 33 OCC mem- 
bers forming the enhancement, and the excess would be ex- 
tended somewhat outward of 0.65 kpc. Below, we show that, 
indeed, we found another kinematic group, the Perseus-Auriga 
group (see Fig.[S]and Table[7J, among the T3 clusters. The com- 
plex is very compact in space as well as in the VPD and consists 
of 8 clusters including one cluster at d xy - 0.77 kpc i.e., out- 
side the main density excess. Consequently, the final census of 
the T3 sample within dxy = 0.65 kpc is 27 OCC 2 members, 
7 members of the Perseus-Auriga group and 19 field clusters. 
This result coincides well with the number of 20 field clus- 
ters estimated before. Although OCC 2 and the Perseus-Auriga 
group are of the same ages, they show different kinematics and 
spatial distribution. Therefore we consider these two groups as 
separate entities. 



lable 4. Candidates of the OCC 2 

COCD Name X Y Z log/ ¥, 

kpc kpc kpc yrs % 



The average age of lcr members of OCC 2 is log t - 8.48, 
with a standard deviation o"i ogf = 0.09. The distribution of the 
T3 clusters in 3D-space is shown in Fig.[8]in the Galactic carte- 
sian coordinate system XYZ. Kinematic members of OCC 2 
form a ring-like structure with an outer diameter of about 1.3 
kpc, and there are only a few T3 clusters in its inner part. The 
ring is ripped along the main Galactic meridian, or more ex- 
actly along a line at I as 15°, / * 195°, where no lcr members 
are found. It even seems that lcr members form two subgroups 
which are separated by this line. For convenience, we call these 
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Fig. 8. Spatial distribution of T3 clusters in 3D-space. The coordinates are given in the Galactic cartesian system. The larger 
circles indicate possible members of OCC2, the diamonds mark members of the Perseus-Auriga group (in order to avoid an 
overcrowding in the X, Z and Y, Z panels, they are shown together with the corresponding boundary circle only in the X, Y panel). 
The dashed circle gives the inner edge of the enhancement area. The other designations are the same as in Fig.0 



subgroups OCC 2a (195° < / < 15°) and OCC 2b (15° < / < 
195°). Comparing the visible distribution of the OCC 2 mem- 
bers with the extinction maps, one can find an agreement in the 
locations of the empty region and the Ophiuchus and Perseus- 
Auriga clouds. Thus, we cannot exclude that some OCC 2 
members are possibly hidden by obscuring clouds. 

Although less prominent, OCC 2a and OCC 2b differ also 
in their location in the YZ plane. OCC 2b shows a dispersed 
distribution and tends to a mean location below the Galactic 
plane, whereas the majority of the OCC 2a members are above 
the Galactic plane. In the VPD of tangential velocities, the 
OCC 2b members are distributed more compactly around the 
centre than the members of the OCC 2a subgroup. 



In Table|4]we list all OCC 2 candidates i.e., all T3 clusters 
located within d n , = 0.65 kpc, which are not Icr members of the 
Perseus-Auriga group. For each cluster, the table contains the 
COCD number, the cluster designation, (X, Y, Z) coordinates, 
the cluster age, and membership probability based on the dis- 
tribution of tangential velocities in the VPD. The OCC 2b can- 
didates are marked by asterisks. 

For eight Icr members of OCC 2 radial velocities are 
known, so we can determine the space velocity vector though, 
less reliably than for OCC 1 (see Tabled). According to the ro- 
tation velocity around the Galactic centre, OCC 2 outruns the 
open cluster centroid with a velocity V - 3.4 +1.2 km/s. The 
average tangential velocity in longitude T/ = 4.0 ± 0.8 km/s in- 



lable 5. Candidates or the Perseus-Auriga group or open clus- 
ters 
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dicates a significant rotation of OCC 2 with respect to the cen- 
troid of open clusters. With an average heliocentric distance of 
OCC 2 of about 0.5 kpc, we find an angular rotation velocity 
a) — 8 + 1.6 km/s/kpc. 

Summarising the main attributes of OCC 2, we conclude 
that the OCC 2 members show a number of properties expected 
for a real physical system. The complex is limited in space and 
it shows a significant (by a factor of 2) density enhancement 
over the surrounding background. The T3 sample has struc- 
tural parameters (Zq and hz in Table Q which differ from the 
other samples of clusters, and for the OCC 2 members these 
parameters are even different from that of the T3 field clusters. 
The OCC 2 members have similar ages and they show simi- 
lar (and relatively compact) kinematics. Therefore, we believe 
that OCC 2 with its 27 lcr-members represents a real group of 
clusters formed about 300 Myr ago. 

The other group of T3 clusters, which we call the Perseus- 
Auriga group due to its special location in space is clearly sepa- 
rated in Figs. [6] and [S] The group has very compact kinematics, 
the average motion being (T/, Tt,) = (-12.6 ± 0.8,5.9 + 1.2) 
km/s with standard deviation of (2.4 + 0.6, 3.5 ± 0.8)km/s. This 
velocity dispersion is even smaller than the cosmic dispersion 
predicted for clusters of corresponding ages by eq. (II Oi and 
shown in Fig.|6] However, the compactness of the velocity dis- 
tribution is not surprising for the Perseus-Auriga group. The 
velocity ellipses derived by eq. dlOi define a measure of cos- 
mic dispersion for clusters of a given age, provided that the 
clusters are distributed all over the sky. Due to the definition of 
tangential velocities, their dispersion depends on galactic lon- 
gitude. For a sample of clusters located within a small range 
of galactic longitudes (like the Perseus-Auriga group), the Tf 
distribution should be more compact than predicted by eq. dlOi . 

Despite of attempts to find more members of the Perseus- 
Auriga group outside the density excess, we could add only 
one more cluster. Thus, the group remains very compact in 3D- 
space. However, this does not exclude the possibility that there 
could be more members since the group is located behind the 
Perseus-Auriga complex of molecular clouds, and we see the 
its members in transparency windows within the clouds. 

Only for two clusters of the Perseus-Auriga group radial ve- 
locities are available. Therefore, the accuracy of the space ve- 
locity is very poor. The data on the Perseus-Auriga group is pre- 
sented in Tableland Table[7] Since all clusters of the Perseus- 
Auriga group have tangential velocities distributed within the 



lable 6. Candidates of the Hyades moving group of open clus- 
ters 
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corresponding dispersion ellipsoid, their kinematic probabil- 
ity is 100%. Therefore, we exclude the probability column in 
Tabled 

5.4. The Hyades moving group of open clusters 

According to Fig.^ the T4 sample does not show any peculiar- 
ity which would require a dedicated discussion. However, this 
sample, also, includes a cluster group which becomes evident 
due to its prominent motion detected in §|4](see Fig. [5}- 

We found six open clusters with considerably large tangen- 
tial velocities. Only one of them, NGC 6991, has a radial ve- 
locity measured, so for this cluster we can determine directly 
the complete vector of space velocity corrected for the Galactic 
rotation which turns out to be (U, V, W) = (-58,-4,2) km/s. The 
other clusters are located at longitudes from 44° to 124° and 
therefore, their —Ti velocities should more or less well coincide 
with U component, which we have estimated to be between - 
30 and -50 km/s. All these clusters are relatively old, with log t 
between 8.57 and 9.16. 

In order to explain a possible nature of their large tangential 
velocities, we considered distribution of the tangential veloci- 
ties corrected for the Galactic rotation and for the Solar mo- 
tion with respect to the cluster centroid versus longitude. Since 
these clusters have large space velocities, there is a good chance 
to observe a sine dependence described by eq. dl 11 1 if the clus- 
ters belong to a common stream. Fig.|9]shows this dependences 
for the T4 sample (panel c). For comparison, the distribution is 
also given for the samples Tl (panel a) and T3 (panel b) where 
the members of the complexes OCC 1, OCC 2 and the Perseus- 
Auriga group are marked, too. 

The visible distribution of tangential velocities of these 
clusters can be fitted by a sine function. This indicates a pos- 
sible stream of several clusters moving in the direction of the 
Galactic anticentre with a velocity of about 40 km/s with re- 
spect to the general cluster centroid. We found three other old 
clusters with large {/-motions which fit to the sine function 
quite well though their tangential velocities are small (due to 
their close location to the centre - anticentre direction). Two 
of them are well known clusters, the Hyades and Praesepe, 
with an age of logf = 8.9 (both) and {/-components relative 
to the cluster centroid of -34.1 and -30.9 km/s, respectively. It 



Fig. 9. Tangential velocities of clusters versus Galactic longi- 
tude. The velocities are corrected for differential Galactic rota- 
tion and for Solar motion with respect to the cluster centroid. 
The points mark clusters of the corresponding sample. Panel a: 
the Tl sample, circles are for the OCC 1 members. Panel b: the 
T3 sample, circles and triangles indicate OCC 2 and Perseus- 
Auriga members, respectively. Panel c: the T4 sample with can- 
didates (circles) of a moving group. The Hyades and Praesepe 
are marked by squares. The distribution of tangential velocities 
of candidates is fitted by a sine function. 

is possible (but not necessary) that these nine clusters form a 
common kinematic stream. The peculiarity of the spatial dis- 
tribution of the clusters prevents us from applying the member 
selection techniques used for other groups. For a clear proof ra- 
dial velocities would be essential. Therefore, we introduce this 
moving group only as a working hypothesis to explain the inho- 
mogeneity in peculiar motion of clusters and call it the Hyades 
moving group after the most famous potential member. This 
cluster group must not be confused with the Hyades superclus- 
ter proposed by Eggen ( 119581 11975b and recently studied by 
Famaey et al. ( 2005 1 who considered several hundreds of stars 
in the Solar vicinity (up to 500 pc) which have similar space 
velocities as the Hyades. In our case, there are a few old open 
clusters which are located at distances up to about 1.2 kpc from 
the Sun but have ages and f/-velocities comparable with those 
of the Hyades. 

The list of candidates of the Hyades moving group is given 
in Tableland Table|7] 

5.5. Summarising cluster complexes 

We classify the complexes OCC 1 and OCC 2 as real phys- 
ical formations, the identification of which is supported by a 
sufficient amount of data. However, it is possible that not all 
their members are completely revealed. For the other two (the 
Perseus-Auriga and specially, the Hyades group), only frag- 
mentary data are available. These make the discussion and ar- 
gumentation of their nature rather difficult. To make a differ- 
ence in the identification status, we refer to them as the groups 
of clusters. In Table0we summarise the parameters of the clus- 
ter complexes and the groups of clusters found in this study. 
The first column gives the names introduced in this paper, N is 
the number of probable members. The average age with the 



Fig. 10. The relative position of identified open cluster com- 
plexes and groups in the XY plane on a background of interstel- 
lar clouds revealed from data on our sample of clusters (Fig.|2j. 
Large symbols mark candidate members for the Gould's Belt 
complex (triangles), OCC 2 (circles), the Perseus-Auriga group 
(diamonds), and the Hyades group (crosses). Small open circles 
mark field clusters. Large circles and the ellipse are apparent 
complex boundaries. 



standard deviation is given in column 3. For OCC 1 which 
probably shows a real age spread, we provide the age span 
only. The next three columns give the spatial data, the coor- 
dinates of the centre (X, Y) and the size of the area occupied 
by the members. In this column, we provide the semiaxes for 
OCC 1 as described above. For OCC 2 the two numbers stand 
for the radii of the outer and the inner borders of the area of 
density excess. For the Perseus-Auriga and Hyades groups, the 
sizes are the radii of areas containing the selected candidates. 
Ti and T/, are the mean tangential velocities corrected for dif- 
ferential Galactic rotation and Solar motion with respect to the 
cluster centroid. For OCC 1, OCC 2, and the Perseus- Auriga 
group they were used for the kinematic selection of members 
(see § 15. 1> . For the Hyades group we computed T/ and Ty only 
with clusters closejo / =_90° (/ = 44 - 124°) where T, * -U. 
The columns U, V, and W give the components of space ve- 
locity corrected for differential Galactic rotation and Solar mo- 
tion with respect to the cluster centroid, and n is the number of 
clusters used for computing these velocity components. This 
number indicates, that only for OCC 1 the space velocity was 
computed with all identified complex members. In OCC 2, less 
than 30% of its members can contribute to the determination 
of the space velocity. For the other groups, these values are 
merely illustrative. Fig^| shows relative location of the dis- 
cussed cluster groups in the XY plane. 

Due to a relative proximity to the Sun and due to the com- 
pleteness of kinematic data, OCC 1 is the most prominent struc- 
ture found in our sample. One of the important features of this 



complex is its clustering in 3D-space. OCC 1 represents a den- 
sity enhancement which is about a factor three higher than the 
surrounding density of field Tl clusters (75 vs. 25 clusters per 
square kpc). It is clearly outlined in space and shows a shape 
which is usually associated with Gould's Belt. It is worth to 
note that up to now Gould's Belt is studied only in connection 
with field OB stars (Frogel & Stothers lT977l Torra et al. l2000l . 
local gas (see Poppel 119971 for a detailed review), and local 
associations (Stothers & Frogel 119741 de Zeeuw Q999 and 
Bobylev 2004 1, and there was no attempt to study it in relation 
with young open clusters in a systematic way. 

Based on the present data we can regard the Gould's Belt 
complex of open clusters as a physical system containing more 
than 20 clusters, occupying a flat elongated area of size of about 
0.8x1x0.2 kpc. The cluster ages range from a few to about 80 
Myr. About half of the clusters are concentrated in a promi- 
nent group residing in the region of the Orion association. They 
have the youngest ages (less the 25 Myr) among the OCC 1 
clusters. The other clusters are distributed along the complex 
more uniformly, and as a rule, they are older than 25 Myr. From 
their location, we find them at the periphery of local molecular 
clouds with embedded associations which expand from Vela 
to Perseus via Scorpio-Centaurus. Only at the top of the spa- 
tial distribution above the Galactic plane (in the region of the 
IC 4665 and Stephenson 1 clusters), OCC 1 emerges out of 
the molecular clouds. The complex shows kinematics typical 
for young clusters, though it is more compact in the VPD than 
field clusters. 

The first hint on possible cluster complexes was made by 
Lynga ( 1982 >, who identified three areas of high concentra- 
tion of young (logf < 7.5) clusters in the Perseus, Carina, 
and Sagittarius regions. Based on comparative statistics of pop- 
ulations of the suggested complexes and of the field, Janes 
et al. (1988 1 found that the three complexes contain mainly 
clusters younger than 40 Myr, and cover areas of order of 
1.2... 1.6 kpc. An extended search for clustering was under- 
taken by Efremov & Sitnik (1988 1. Among stellar associa- 
tions and young clusters (earlier than b2 in the scale of Becker 
& Fenkart I197U . they identify 17 groups of various sizes 
(150... 700 pc) containing at least 3 objects. They suggested that 
complexes of young clusters both outline the spiral arms and fill 
inter-arm space in the Solar neighbourhood (see Efremov 1995 
for more details). 

The presence of two cluster complexes and two more clus- 
ter groups in the Solar neighbourhood, with only one being re- 
garded as a young object implies that such aggregates seem to 
be quite usual in the galactic disk. Therefore, it is not necessary 
to relate them to spiral arms or current star formation. If we do 
not suggest an exclusive role for the Solar neighbourhood, we 
may expect to find footprints of cluster groups (or their rem- 
nants) of different age everywhere in the Galactic disk. With 
time, such groups should loose their members due to cluster 
evaporation, become less prominent and will be detectable in 
the velocity space only. 

We would like to stress that revealing the complexes and 
groups of clusters was not an output of a systematic search for 
such objects in the Galactic disk but rather the result of an at- 
tempt to understand the nature of irregularities either in density 



or in kinematic distribution which we found within the com- 
pleteness area of the cluster sample. Therefore, one may expect 
that a systematic search in larger samples (which perhaps are 
less complete but cover larger areas) will result in identification 
of additional cluster agglomerations. 

6. Formation rate and lifetime of clusters 

One of the first extensive studies into the distribution of ages 
of open clusters and the estimation of their lifetime was car- 
ried out by Wielen dl971> . The statistics of cluster ages was 
based on data from the catalogues of Lindoff |l968l, and 
Becker & Fenkart (I197U . Later studies on this topic by Janes & 
Adler fE?52l . Lynga ( H9521 . Pandey & Mahra ( fT986l . Pandey 
et al. fi987l . Janes et al. JT9881 . Bhatt et al. (1989} and 
Battinelli & Capuzzo-Dolcetta ( 1991 1 used the approach pro- 
posed by Wielen ( I197U . but considered the larger list of clus- 
ters by Janes & Adler ( 1982 1 as well as the third, forth, and fifth 
versions of the Lund catalogue (Lvnga ll983l IT9851 \TW7\ . In 
comparison to these studies, our sample contains at least twice 
as many clusters with ages and distances determined by a uni- 
form technique, and it includes many newly detected clusters. 
Therefore, a revision of the previous findings and conclusions 
is reasonable and relevant. 

Let dN(t) = if/(t) dt be the number of clusters formed in 
some interval of time (/, / + dt), where iff(t) is the cluster forma- 
tion rate (CFR). Further, we assume that the clusters decay with 
a probability p = 1 /t where r is the typical lifetime of these 
clusters. If N(a) is the number of clusters with age a, which are 
observed at the present time t\ — t + a, we can relate N(a) to 
the CFR via the equation 

N(a) = N(h - a) e~ fl/T = if/fa - a) e~ fl/T . (13) 

If clusters would never dissolve (i.e., r = oo), then their age 
distribution would reflect the temporal behaviour of the cluster 
formation rate. In contrast, for a constant CFR the age distri- 
bution N(a) is controlled by the decay process. In the general 
case, the observed distribution of cluster ages depends on both 
factors, but in principle, it is possible to fix the processes. Since 
the observed span of cluster ages (a few Gyr) is much shorter 
than the age of the Galactic disk (over 10 Gyr) and since no 
indications of strong temporal variations are observed for the 
local rate of star formation (see e.g. Scalo 1986 1, it is rea- 
sonable to assume that the formation rate of clusters is time- 
independent. Analysing the age distribution of open clusters 
under the assumption of a constant CFR during the disk life- 
time, Wielen (119711 found a cluster half-lifetime t\/2 of 160 
Myr, and a cluster formation rate of - 0.1 kpc~ 2 Myr~ 1 . Here 
fi/2 is the decay time of one half of the clusters once formed, 
and it can be derived from eq. J 1 31 as t\/2 — — t In 1/2 « 0.7 r. 
Using larger samples of clusters, the authors above obtained 
about the same or somewhat smaller lifetimes of open clusters. 
From a comparison of the age parameters for clusters of dif- 
ferent richness classes and of different spatial location in the 
Galaxy, a reasonable conclusion had been achieved that the 
cluster lifetime increases with the richness class as well as out- 
wards from the Galactic centre. 



lable 7. Parameters or newly identified cluster complexes and groups 



Complex, N \ogt X c Y c size T, T b U V W n 

group yr kpc kpc kpc km/s km/s km/s km/s km/s 

OCC 1 23 6.67-7.88 -0.078 -0.053 0.385,0.475 -3.6 ± 1.4 +1.0 ± 1.0 -4.0 ± 1.2 -0.2 ± 1.1 +0.3 ± 0.6 23 

OCC2 27 8.48;0.09 0.000 0.000 0.325,0.65 +4.0±0.8 +2.9 ± 0.7 +1.9±4.2 +3.3 ± 1.2 +2.5 + 1.3 8 

Perseus-Auriga 8 8.47;0.09 -0.450 0.300 0.26 -12.6 + 0.8 +5.9 + 1.2 +4.6±13.6 +21.1 + 9.3 +4.3 + 4.3 2 

Hyades" 9 8.83; 0.21 0.080 0.540 0.66 +34.9 + 3.7 -0.4 + 4.5 -34.4 + 5.3 -0.8 + 3.8 +2.0 + 2.6 4 

1 The mean tangential velocities were computed with 6 clusters (without the Hyades, Praesepe, and ASCC 99) 



Table 8. Cluster formation rates (CFR) and cluster lifetimes 
(t) 



No 


CFR 

kpc~ 2 Myr~' 


T 

Myr 


N d 


Sample 


Note 


1 


0.10+0.01 


256+12 


652 


tot 




2 


0.31+0.03 


327+25 


259 


cmp 




3 


0.23+0.03 


322+31 


195 


fid 




4 


0.10 


231 


70 


w71 


1 


5 


0.18 


144 


112 


p86 


2 


6 


0.25 


100 


213 


j88 


3 


7 


0.45 


14 


100 


b91 


4 



Sample: No 1...3 are the present study, No 4.. .7 are taken from the lit- 
erature; tot - the total sample; cmp - the sample of clusters within the 
completeness area; fid - the sample of field clusters within the com- 
pleteness area; w71 - Wielen 119711 : p86 - Pandey & Mahra 119861 : 
j88 - Janes et al. (19881 ; b91 - Bat tinelli & Capuzzo-Dolcetta (19911 . 
Notes: 1 - the Becker & Fenkart <1971t sample; 2 - the Lund cata- 
logue, 3rd edition; 3 - the Lund catalogue, 5th edition; 4 - clusters 
with d < 2kpc & V,„, fg , < -4.5 mag. 



Fig. 11. Distribution of open clusters versus age. For an eas- 
ier comparison, the distributions for the different samples are 
not normalised to unit area. Upper panel: our total sample 
is shown as the hatched histogram. The sample of field clus- 
ters within the completeness area is marked as the filled his- 
togram, and the solid curve marks the fitted age distribution. 
The age step is 50 Myr. Lower panel: the same distributions 
as in the upper panel together with the age sample used by 
Wielen (I197U . The data from Wielen J1971I are shown as the 
backhatched histogram in the foreground. The vertical bars cor- 
respond to Poissonian errors derived from cluster counts. The 
binning is chosen same as in Wielen dl971> . 

Taking into account the statistical properties of our cluster 
sample, we computed the CFR and the cluster lifetime for three 
cases. At first, we consider the total sample of 652 clusters. In 
the second case, we include only those 259 clusters which are 
located within the completeness area (d xy < 0.85 kpc). Finally, 
in order to eliminate possible biases in the age distribution due 
to the cluster complexes, we excluded their lcr-members from 
the second sample. This sample of so-called field clusters con- 
tains 195 clusters. 

In the upper panel of Fig. the distributions of clusters 
versus age are shown for the total sample as well as for the 
field clusters within the completeness area. As expected, the 



total sample is biased towards young clusters. Compared to the 
distribution of field clusters, this results in a steeper average 
slope of the relation and, consequently, in smaller lifetimes. 
The curve obtained from fitting eq. J 1 3I > coincides well with 
the observed distribution of field clusters. The parameters of 
the fit are given for all three cases in the first three rows of 
Table [8] We conclude that the completeness issue has a strong 
impact onto the determination of CFRs and cluster lifetimes. 
The exclusion of the complex members does not change the 
lifetime t significantly, whereas the impact on i// is somewhat 
stronger: it decreases the CFR by about 20%, but the difference 
is still within the mean errors. 

In the bottom panel of Fig. ^2 we show the same distri- 
butions together with results from Wielen dl971> based on the 
Becker & Fenkart ( 1197U sample. In its shape, the distribution 
used by Wielen dl971> is in better agreement with the data of 
the total sample, but even so, there is a pronounced deficiency 
of older clusters which play an important role in the determina- 
tion of cluster lifetime. This is the reason for the smaller life- 
time derived by Wielen dl97U . given in Table|8] too. 

Table [8] also includes the corresponding findings by other 
authors from the 80s and 90s. To ease comparison, we trans- 
formed 1 1/2 to t, if the cluster lifetimes were not explicitly given 
in the original publication. Not only do the determinations of 
t differ strongly among themselves, they are far off Wielen's 
results, and even farther from ours. We explain the disagree- 



ment or our results with the results or Fandey& Mahra ( 1986 1, 
and Janes et al. ©88), based on the Lund Catalogue, by the 
considerable increase of data on clusters older than 250 Myr in 
the last decade. For example, the number of clusters with ages 
of 250 . . . 3000 Myr within 0.85 kpc from the Sun has been in- 
creased by a factor of 6 due to the newly discovered clusters. In 
comparison, the total number of clusters has grown from 1 150 
in the 5th edition of the Lund Catalogue to about 1700 clusters 
in Dias et al. (2004J i.e., by only a factor of 1.5. Concerning 
the extremely low lifetime t derived by Battinelli & Capuzzo- 
Dolcetta (1991 the selection effect is even stronger. Since they 
considered only the brightest clusters, their sample is overabun- 
dant on young objects, and the lifetime is strongly underesti- 
mated. 

The footprints of selection effects become also visible in 
the results for the CFR in Table [8] Unlike the lifetime deter- 
mination which mainly depends on a realistic relation between 
numbers of young and old cluster, the CFR results are more 
influenced by overall incompleteness of the sample studied. 
The increase of the CFRs computed with our samples corre- 
lates clearly with an increasing completeness of the data (cf. 
Table[S}. A similar correlation can be observed for the CFR re- 
sults derived by Wielen ( H971L Pandey & Mahra dl986> . and 
Janes et al. ( 1988 ). The good agreement between our CFR and 
that from Janes et al. (0.988 1 is expected, since we use the same 
basic source, the 5th edition of the Lund Catalogue. On the 
other hand, we believe that the considerable disagreement in 
the CFR derived by Battinelli & Capuzzo-Dolcetta ( 1991 ) and 
all others is a consequence of the underestimation of ages of 
their young clusters. It is well known that the turn-off age cal- 
ibrations are biased in the case of young clusters (a few tens 
of Myr or younger) due to the rather steep Main Sequences 
of early-type stars towards younger ages. Neglecting this bias, 
one underestimates the age spread within a sample of young 
clusters, and hence, overestimate their abundance. 

Taking the cluster lifetime derived above and an age t\ for 
the Galactic disk of about 10-12 Gyr, we compute the total 
number of cluster generations n g - t\jr as being of the or- 
der of 30-40. Provided that a typical open cluster consists of 
n* = 1000 stars at its birth and with the surface density of clus- 
ters being "L — 114 kpc -2 (see § I3.3> . the total surface density 
n* of disk stars and their remnants which have passed an open 
cluster member phase is n g n* a E « 4 x 10 6 kpc -2 . This number 
is in a good agreement with another estimate coming from the 
derived CFR i.e., n* = ri^iffti » 3.6 x 10 6 . Compared with the 
local density of disk stars of about 7 x 10 7 kpc -2 , we find that 
less than 6% of the total stellar population of the Galactic disk 
is genetically connected with classical open clusters. On the 
other hand, this estimate is based on assumptions and extrapo- 
lations which at present can not be checked in detail. Therefore, 
we can only conclude that that portion of the disk stellar popu- 
lation, which had been originated in open clusters, is relatively 
small, of the order of 10%. 

This estimate is implicitly supported by another statistics 
available from the ASCC-2.5. Let us count all stars which are 
brighter than V — 9 mag i.e., a magnitude range where 99% 
of stars have spectral classification and, practically, where all 
clusters are identified in the ASCC-2.5. Among the intrinsi- 



cally bright population (supergiants, O - B3 stars), we then 
computed the ratio of lcr and 2cr cluster members to field 
stars. It turns out that about 10% of these 3770 extremely lu- 
minous stars are members of open clusters. Our results sup- 
port the "pessimistic" estimate of 10% obtained by Miller & 
Scalo ( 1978 1 from a much poorer statistics. 

The rather low input of open clusters to the stellar popula- 
tion of the Galactic disk can be explained by their high infant 
mortality, reducing generations of newly born clusters by a half 
(Lamers et al. 2005 1. Even so, this is not sufficient to account 
for all Population I stars. Thus, the current statistics support the 
early conclusion by Miller & Scalo ( 1978 1 that about 65% of 
the disk stars should have their origins in the association reser- 
voir. 

7. Conclusions 

In this paper we studied unbiased properties of the local clus- 
ter population taking advantage of a homogeneous set of open 
cluster parameters derived from the all-sky census of cluster 
stars up to V — 12.5. The revised cluster membership, the im- 
plementation of accurate kinematic and photometric data sup- 
plemented by radial velocities and spectral classification avail- 
able at the moment, allowed us to extend the completeness area 
of open cluster data in the optical up to a distance of 0.85 kpc. 
The resulting sample was the basis for the determination of spa- 
tial, kinematical and evolutionary parameters of cluster popu- 
lation in the local Galactic disc. 

Although the spatial distribution of clusters on scales com- 
parable to inter-arm distances (d xv w 2 kpc) is compatible with 
a model of a constant surface density, the distribution is not 
uniform. There are two main factors influencing the observed 
distribution in the optical. At first, the cluster can form groups 
which are revealed in the age distribution as well as in the space 
and velocity distributions. We found three cluster complexes 
in the Solar neighbourhood affecting significantly the apparent 
density pattern. The other effect influencing the apparent dis- 
tribution of nearby open clusters in the optical is a patchy in- 
terstellar extinction due to local screening clouds. Also, these 
have a considerable impact on the observed distribution of clus- 
ters at large distances which we can see only through trans- 
parency windows. It is not excluded that some features of the 
local disk (like the Local arm) which have been discussed for 
a long time, are only virtual details which appear due to an 
uneven distribution of the reddening. The vertical distribution 
of the cluster sample turned out to be extremely flat, with a 
scale height of about 50. . . 60 pc, and it does not depend on 
age at least for cluster ages up to 1 Gyr. We found the sym- 
metry plane of open clusters to be Zo = -22 ± 4 pc, and the 
total density of clusters in the symmetry plane is D(Zq) = 1015 
kpc -3 . The spatial distribution of younger clusters is compati- 
ble with a flat grand design of a spiral pattern with a pitch angle 
of p = -5° . . . - 6°. 

The kinematical parameters of the cluster sample are com- 
parable to those of young field stars in the Solar neighbour- 
hood though, the clusters cover distances which are typically 
larger by a factor of two than Hipparcos-based samples of field 
stars. This confirms that, in general, the observed cluster pop- 



ulation is young. A combination or accurate cluster ages and 
kinematical parameters provides a possibility to study the tem- 
poral variation of the cosmic velocity dispersion. We found that 
in average, the dispersion of each velocity component increases 
by a factor of two during a period of 3 Gyr. The local velocity 
field within about 1.5 kpc can be sufficiently well described by 
the linear model of differential rotation. But similar to the den- 
sity distribution, the regular velocity pattern can be broken by 
groups of clusters with significantly different motions. So, we 
found a group of nine older clusters showing a large velocity 
component in the Galactic anti-centre direction. 

Fluctuations in the spatial and velocity distributions can 
be attributed to the existence of three open cluster complexes 
containing up to a few tens of clusters of similar ages. The 
youngest complex, OCC 1 (logf < 7.9), is a structure of about 
1 kpc in size inclined to the Galactic plane by 19°. Highly prob- 
able, it is a signature of Gould's Belt which has generally been 
considered in the literature to be an aggregation of the youngest 
field stars, gas and associations. The most populated complex, 
OCC 2, includes clusters of moderate age (logf * 8.45). The 
clusters of the Perseus-Auriga group have the same age but 
show different and compact kinematics. They are also very 
compact in their location and seen in voids of the Perseus- 
Auriga clouds. Possibly, we see only a part of a larger struc- 
ture covered by clouds. The oldest (log f w 8.85) and the most 
loose group was detected due to a large peculiar motion (i.e., 
velocities corrected for Galactic rotation and Solar motion) of 
about 25-45 km/s in the Galactic anticentre direction. We note 
that the Hyades and Praesepe have about the same age (logf 
= 8.9) and they have an U velocity component of -34.1 and - 
30.9 km/s, respectively. Nevertheless, the real existence of the 
Hyades moving group of open clusters cannot be proven due to 
the lack of radial velocities of the suggested candidates. 

From the parameters of spatial distribution, we estimated 
a total number of 10 5 open clusters currently in the Galactic 
disk. Cluster lifetime and formation rate obtained from the age 
distribution of field clusters within the completeness area are 
found to be 322+31 Myr and 0.23 ± 0.03 kpc _2 Myr -1 , respec- 
tively. Assuming a typical open cluster of the Pleiades type, one 
derives the total surface density of disk stars passed through 
the phase of open cluster members to be about 4 x 10 6 kpc -2 . 
Compared to the local density of disk stars of about 7 x 10 7 
kpc -2 , one obtains that the input of open clusters into the total 
population of the Galactic disk is about 6%. 
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